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1 Introduction 



The threshold network model is a type of finite random graphs. In this paper, we introduce a 
generalized threshold network model. A pair of vertices with random weights is connected by an edge 
when real-valued functions of the pair of weights belong to given Borel sets. We extend several known 
limit theorems for the number of prescribed subgraphs to show that the strong law of large numbers 
can be uniform convergence. We also prove two limit theorems for the local and global clustering 
coefficients. 

^ • Complex networks have been an attractive research topic for a decade. Particularly, many real-world 

graphs are characterized by the small diameter, high clustering (abundance of connected triangles), and 
fat-tail degree distributions. Degree distributions often follow the truncated power law, which is called 
the scale- free property of networks [1,4,15]. Both growing and static network models are capable of 
^ ' generating scale-free networks. 

I Here we are concerned with asymptotic properties of a class of static network models called the thresh- 

\l ' old network model, which is generated on n vertices labeled 1, . . . ,n with independent and identically 

distributed (i.i.d.) random weights Xi, . . . ,X„. We connect a pair of vertices i and j with « 7^ j by an 
edge when Xi + Xj > for a given threshold 6. The threshold network model is a subclass of so called 
hidden variable models and its mean behavior [5,6,8,10,11,18,19] and limit theorems [9,14] have been 
analyzed. 

\ To define a generalization of the threshold network model, let be the d-dimensional Euclidean space. 

We prepare an i.i.d. sequence of R'^-valued random variables Xi, . . . ,Xn with a common distribution 
function F. We associate the random variable Xi, which we call weight function, with vertex i. Now we 
^ . introduce Borel measurable functions : (M'')^ R with fj: (a;, y) = fj, (y, x) for all I' € {I, ... , I}. Let 

H I B(R) be the Borel ct- field of M. For a given finite collection of Borel measurable sets C = {Bi, . . . , Bi} 

■ with Bii e B{R), we connect vertices i and j {i ^ j) if f^^{Xi,Xj) e Bj' for aU I' e {1,...,^}. In 

other words, we form an edge if Y\\,^iIbi, (^fl {Xi,Xj)^ = 1 for i 7^ j, where denotes the 

indicator function, i.e., Ia{x) = 1 for x € A and Ia{x) = otherwise. Thus we obtain a random graph 
Gc{Xi, . . . , Xn). If there exist two collections of Borel sets C = {Bi, . . . , Bi} and C = {B[, . . . , B'l} with 
Bi, C B'l, for all I' G {1, . . . ,/}, then P{(z, j) e Gc{Xi, . . . ,X„)} < P{(i,j> e Gc'(^i, • ■ ■ ,^«)} holds by 
simple coupling. 

This random graph generalizes the threshold network model studied in [5,6,8-10,12,14]. By choosing 
I — 1, Bi = {9, 00) for some G M, /c (a;, y) = x + y, we reproduce the model in [5,6,8-10]. In the context 
of social networks, a model with I — 2, Bi ~ {9, 00), B2 = (— oo,c] {9,c S R), fl[x,y) = x + y, and 
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/^(a;,y) = |x — y| (or f'^{x,y) = \x — y\/{x + y)) was proposed [12]. General limit theorems are shown 
in [14] when I = 1^ Bi = {6, oo), and f^{x, y) = \x — y\. 

In Sec. 2, we state several general limit theorems for the number of prescribed subgraphs. By using 
[/-statistics, the strong law of large numbers, the central limit theorem, and the law of the iterated 
logarithm are stated for global properties of the model. We also state a limit distribution for a local 
property. These are generalizations of Theorems 1, 4, and 5 of [9] and Theorems 1, 2, and 3 of [14]. In 
Sec. 3, we show that the strong law of large numbers for the number of prescribed subgraphs is uniform 
convergence on so-called the VC class of Borel sets, generalizing Theorem la of [14]. In Sec. 4, we show 
limit theorems for the clustering coefficient, which quantifies the abundance of connected triangles in a 
graph in a specific ways. Particularly, we show the strong law of large numbers for the local clustering 
coefficient (Theorem [5]) and the global clustering coefficient (Theorem [3]) . Theorems [2] and [3] are main 
results of this paper. In Sec. 5, we present several examples of limit degree distributions. 



2 General Results 

In this section, we show limit theorems for the number of prescribed subgraphs. Let us begin with 
notations [14]. For m e {2, . . . , n}, we consider a graph H = {Vh, Eh) on the ordered set of m vertices 
Vh = {vi, . ■ ■ ,Vm) and the edge set Eh- For another graph H' ~ {Vh'jEh') on m vertices, we say 
H' ~ H a Vh' = Vh and Eh' — Eh for some reordering of vertices. Thus = {H' : H' ~ H} is the 
set of all graphs isomorphic to H. Let us define Am — Ui -^ff ' where Hi is an arbitrarily chosen graph 
on m vertices. The collection of all triangles and graphs on three vertices that consist of two connected 
vertices and an isolated vertex is an example of ,43. The collection of cliques on m vertices and the 
graphs on m vertices with m isolated vertices is an example of Am- The definition of Am guarantees the 
symmetrical property of the kernel function hj\^ : (M'')™ — > M given by 

ftyt„(a;i, . . -,Xm) = lAr„ {Gc{xi, - - -,Xm)) , (1) 

where Gc{xi, . . . , Xm) denotes a realization of the random graph Gc{Xi, . . . , Xm)- Then we define 

tjn{C,Am)— ^ fT'A,„{Xtj^, - - - , Xi^), 

l<ii <---<ij^i<n 

i.e., the number of subgraphs belonging to the collection Am in the random graph Gc{Xi, . . . ,X„). We 
also define 

Un{C,Am;i)— ^ f^A^iXi, Xi^, . . . ,Xi^), 

I<i2<---<i77i<n 

i.e. , the number of subgraphs that include vertex i and belong to Am in the random graph Gc {Xi , . . . , Xn)- 
Notc that Un{C, Am] *)i 1 — * — identical in distribution and the following relation holds: 

Etl UnjC, Am-, _ mUn{C,Am) _ UnjC , Am) 

<m-\) O ■ 

This implies that the global property Un{C,Am)/{J^ is the arithmetic mean of the local properties 

t/„(C, 1)/CZ\) , • ■ ■ , f/«(C, A„; • 
We define 

F{C,Am)^nhA^{Xl,---,Xm)l 

aC^Am) = Var{nhA„AXu - - - .Xm)\X^]), 

and assume C,{C,Am) > 0. Since Un{C , Am) / (^) is a [/-statistic [17] obtained from the symmetric kernel 
^Amj the strong law of large numbers (SLLN), the central limit theorem (CLT) and the law of the iterated 
logarithm (LIL) are derived from general results for the [/-statistics, namely. Theorem A (SLLN) and 
Theorem B (LIL) in Section 5.4, and Theorem A (CLT) in Section 5.5 of [17]: 
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Fact 1 (SLLN for global property). As n oo 



F{C,Am), almost surely. 



Fact 2 (CLT for global property). As n ^ oo, 



z, 



m2C(C, A 

where =^ stands for convergence in distribution and Z is a standard normal random variable. 
Fact 3 (LIL for global property). 

Un : •Am ) 



' n(log log n)~i 



limsup,. 



(") 

\7n/ 



i^(C,An) 



1, almost surely. 



There are the direct generalization of Theorem 4 of [9] and Theorems 1, 2, and 3 of [14] to the present 
model. 

Remark 1. Generally, when the fj: is asymmetric (e.g. directed graph) , the number of graphs isomorphic 
to a graph H on m vertices is ("j) • ml. The limit theorems above are valid by replacing the normalizing 
factor {::;) wUh{-)-ml. 

By generalizing Theorems 1 and 5 of [9], we obtain the following asymptotic behavior: 

Fact 4. As n ^ oo, 

Un : Arm ^ ) 



Vm-1/ 



U{C,A 



where 



U{C,Am) 



hA^{Xl,X2, . . .,Xm)F{dx2) ■ ■■F{dXra)- 



Remark 2. Almost sure convergence theorem for Un{C,Am] sense of Theorem\^in Sec. 

4-1 can be proved by a simple modification of the proof of Theorems 1 and 5 of [9]. 



3 Uniform Property 

The Vapnik-Chervonenkis approach is well known in the context of the statistical learning theory. Par- 
ticularly, it is useful in showing uniform convergence for limit theorems [7, 16]. In this section, we show 
that SLLN for global property (Fact [1]) is uniform convergence on the VC class of the Borel sets, which 
extends the special case treated in [14]. 

Let M be a set and P be a family of subsets of A-f. For A C M let A'^(A) =^{Ar^V), where ^{Ar^V) 
denotes the number of sets in A n I? = {A Ci D : D £ T>}. Let mx>{n) = max^cM {^'^{A) : \A\ = n} for 
n = 0, 1, 2, . . ., where \A\ denotes the number of elements in A, or if |M| < n let m-p(n) = mx> {\M\). We 
define an indicator of the family T>: 

S{T>) \ ^^P^"" ■ "^f ("■) ~ 2"} if T) is non-empty, 
1—1 if I? is empty. 

The family T> is called a Vapnik-Chervonenkis (VC) class of sets if S{T>) < +00. For example, the 
collection of half intervals T> = {{—oo,x] : a; G R} is a VC class on R with S{T>) = 1. Based on Chapter 
4.5 of [7], we have 
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Corollary 1. For any V C 2^' and V C 2^' (resp. 2^), if S(V) < oo and S{V') < oo then 
S{V \JV') <oo and S{V C^V') < oo (resp. S{V xV) <oo) 

where VUV = {DU D' : D gV, D' G V'}, V nV = {DClD' : D eV, D' eV'} and V xV = 
{D X D' : D e V, £>' e V'}. 

For a given function h : M R, the subgraph of h is the set 

Dh = {(x,t) € M xR:0 <t < h(x) or h{x) <t<{)} . 

A class of functions 7i is a VC-subgraph class if the collection V-h — {Dh : h E H} is a VC class of sets. 

For a class of real-valued measurable functions H on Af™ for a fixed integer m, Arcones and Gine [2] 
proved the following uniform SLLN for i.i.d. sequence {Xi}^^^ ^ on M: 

Lemma 1 (Corollary 3.3 of [2]). //7i is a measurable VC-subgraph class of functions with 
E [sup^jg^ \ h{Xi, . . .,Xm)\] < +0O, then 



sup 

hen 



\rnJ 1 



E 

<«1 < ■ ■■<irn 



h{X,„...,X,J-E[h{X,,...,X^)] 



0, almost surely, 



oo. 



In order to use Lemma [1] we rewrite the kernel function ([T]) and show that the family of the kernel 
functions indexed by the VC class of Borel sets is a VC-subgraph class. We assume d = I = 1. For a 
fixed integer m > 2 and : ^ M, we define a function G : —> rC?) by 

G{xi, . . .,Xjn) = (/c(a;i,X2),/c(a;i,a;3), . . . , f}{xi,x„i), fXx2,X3), . . . , f}{xjn-i,Xm)) ■ 

Each coordinate corresponds to a potential edge of the graph with a lexicographic order. For example, 
if m = 4, the first coordinate corresponds to the (1,2), the second coordinate to (1,3), and the sixth 
coordinate to (3,4). Note that edge {i,j) exists if and only if f^{xi,Xj) S B for a given Borel set B. 

For a given collection Am of graphs on m vertices, we define a set Am on k( ^ ) as follows. For each 
graph H G Am, we associate a set -ff on r( ^ ) _ When a pair of vertices in H has an edge, the corresponding 
coordinate of H is occupied by B. Otherwise, it is occupied by B'^, where B'^ denotes the complement 
of B. For example, if m = 4 and has edge set {(1, 2), (2, 3)}, then H = B x B" x B" x B x B" x B" . 
Then we define the set Am = V}h£A ^- Note that G{xi, . . . , Xm) G Am is equivalent to the event that 
the realization of the random graph with weights xi, . . . ,Xm is in Am- Finally, we obtain the rewritten 
form of the kernel function: 



[xi 



(G'(xi, . . .,Xm)) 



.(xi,.. 



0, 



(2) 



where G^^{Am) denotes the inverse image of Am- 

Let I? be a VC class of Borel sets on M. For a fixed collection Am, we consider the class of kernel 
functions Ti-p = {/i^^(a;i, . . . , Xm) '■ B G I?}. By Corollary [U if I? is a VC class of sets in general, then 
the class of indicators {Ijj : D e V} is a VC-subgraph class. From Eq. ([2]), if the collection G^^{Am) = 
{G~^{Am) ■ B E V} is a VC class then T-C-d, which is a collection of indicator functions, is a VC-subgraph 
class on R™. Indeed, G^^{Am) is a VC class based on Theorem 4.2.3 of [7] and Corollary [T] above, and 
Hv is a VC-subgraph class. Finally, we find the uniform version of Fact [T] by Lemma [H 

Theorem 1. IfD be a VC class of Borel sets, then for a fixed fc and Am, 



sup 



Vm/ l<ji <■■■<»,„<« 



,X,J~¥.[h^JXi,...,Xm)\ 



almost surely, as n 



Theorem[T]can be extended to general d and I. It generalizes the uniform SLLN in Theorem la of [14], 
which deals with the collection of half intervals as a VC class of sets. 
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4 Clustering Coefficient 



Real-world networks are often equipped with high clustering, that is, a large number of connected trian- 
gles. The clustering coefBcicnt quantifies the density of triangles in a graph (see [1,15] for review). In 
this section, we study the limit theorems for the clustering coefficient. 



4.1 Local Clustering Coefficient 

We assume d = I = 1; extensions of the following results to general d and / is straightforward. We 
consider a random graph Gb{Xi, . . . , X„) for a given Borel set B and fc = fl- Then we define 

l<j<n 

Tn{i) = ^ hT{Xi,Xj,Xh), 

l<j<k<n, 

where hD{x,y) = lB{fc{x,y)) and hT{x,y, z) = lB{fc{x,y)) ■ IsUdy^z)) ■ lB{fciz,x)), i.e., D„(i) is the 
degree of vertex i and T„(z) is the number of triangles including vertex i. The local clustering coeflBcient 
C„(z) of vertex i is given by 

'^"(*) = "^lyy ■ hDn{i)>2} + w ■ /{D„(i)=o,i} 

for an indeterminate w. The second term represents the singular part for which the local clustering 
coefficient is not defined in physics literature and applications. Here we retain this term to assess the 
contribution of vertices with degree or 1. If it is necessary to restrict C„(i) G [0, 1], we must substitute 
a real value on [0,1] into w. If we substitute into w, the contribution of these vertices to C„(i) is 
ignored. If we substitute 1, this contribution is implied to be the maximum because vertices with degree 
more than one satisfies C„(i) < 1. Now we define 

Vn{i)= hv{Xi,Xj,Xk), 

l<j<k<n, 

where hv{x,y,z) = lB{fc{x,y)) ■ lB{fc{x, z)), which represents the number of vertex pairs (j, fc) such 
that both vertex j and vertex k are connected to vertex i. We note the relation: On {Dn{i) > 2} or 
equivalently {Vn{i) > 1}, 



Dn{i) 



Vnii), 



which leads to 



We also define 



where 



Tnii) 



Et {X) 

= En{x')'^ ^{Ed{x,)>o} + w ■ I{Eo{x,)=o}, (3) 



ED{Xi)= [ hD{Xi,y)F{dy), 
ET{Xi)= [ [ hT{Xuy,z)F{dy)F{dz). 

JR JR 



Wc consider Cn{i;x), C{i;x), Dn{i;x), Tn{i;x), and Vn{i;x), which are random variables Cn{i), C{i), 
Dn{i), Tn{i) and Vn{i) restricted to the subspace such that {Xi = x}. For example, T„(l;a;) = 
^2<j<k<n^T{x,Xj,Xk). We obtain the following asymptotic results for C„(i): 
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Theorem 2. As n oo, 

(i) For any a: S M, C„(l;a;) — ^ C{\;x), almost surely. 
In particular, 

(ii) C„(l) —^ C{1), almost surely. 

Proof. For an arbitrary fixed a; e R, we first prove 

E[hvix,X2,X3)] = ^P{Vnil;x) = for all n > 1) = 1. (4) 
Indeed, if E [V(a;, ^2, ^3)] = then 



[Vn{l;x)]=: 



hv{x,Xj,Xk 

2<j<k<n 



= nhv{x,X,,Xk)]=0 

2<j<k<n 



for all n > 1. Conversely, if E [14 (1; x)] = for all n > 1 then 

E[hv{x,X2,X3)]=E[V3{l;x)]=0. 

Therefore, we obtain 

E[hv{x,X2,X3)] = <==>E[Vnil;x)] = for aU n > 1. 

Since Ki(l;a;) is nonncgative, 

E[T4(l;a;)] = P (F„(l; a;) = 0) = 1 for aU n > 1. 

Moreover, {V^(l;a;) — 0} is nonincreasing with n, which implies 

F{Vnil;x) = 0) = 1 for aU n > 1 P(F„(l;x) = for all n > 1) 1. 

Thus we have Eq. (U. By definition, Vn{l;x) is invariant under any permutation on {a;2, ^3, . . . , a;„}, 
and Vn{l',x) is nondecreasing, i.e., 

K(l;a;)(a;2,a;3, ...,Xn)< T4+i(l; a;)(a;2, ^3, . . . ,Xn,Xn+i) (5) 

for all n > 1. Therefore P x) = for all n > 1) equals to zero or one by the Hewitt-Savage zero-one 

law (see Theorem 36.5 of [3]). So we have 

E[hv{x,X2,X3)] > P(V;i(l;a;) = for aU n > 1) P (Kill; > 1 for some n > 1) = 1. 

Using Eq. (O, {Vn{l;x) > 1 for some n > 1} is equivalent to the event 
{3N > 1 s.t. Vn{l;x) > 1 for all n > N}. Hence we obtain 

E[hvix,X2,X3)] > F{3N > 1 s.t. Vnil:,x) > 1 for aU n > A^) = 1 

^ P{3N > 1 s.t. C„(l;a;) = T„(l; for all n> N) ^1. (6) 

Since hxix, X2t X3) and hv{x, X2, X3) are symmetric functions of a;2 and X3, we define i7-statistics 

/i-'n ^ ^ 7;r^ hT{x,Xj,Xk), 

\ 2 ) \ 2 ) 2<j<k<n 

— 7— — TT— - — TT~ ^ hv (x. X^ ^ Xk). 

12/ \ 2 I 2<j<k<n 

We have the following SLLN by Theorem A in Section 5.4 of [17]: As n — > 00, 

Tn{l;x) 



Vn{l;x) 



E[hT{x, X2, X3)] , almost surely, (7) 



E[hv{x, X2, X3)] , almost surely. (8) 



Based on Eqs. ([7]) and the corresponding clustering coefficient 

r„(i;x) _ rn(i;x)/("-^) 

converges to E [hT{x, X2, X^)] /E [hv{x, X2, X3)], almost surely as n ^ 00. By Eq. we have 

lim Cn{l;x) =E[hTix,X2,X3)]/E[hv{x,X2,X3)]) ^ I. (9) 



On the other hand, Eqs. (HI and © imply that E [hv{x, X2, X3)] = is equivalent to 

P (lim„^oo Cn{l; x) = w) = 1. With the relation E [hv{x, X2, X^)] = E [hi:i{x, ^2)]^, we obtain 

r n \ r^^ \ nhT{x,X2,X^)] 

E[hD{x,X2)\ 

almost surely as n ^ 00. Particularly, we have by using Fubini's theorem, 

Pf lim C„(l) =C(1)) = [ P( lim Cn{l;x) = C{l:,x))F{dx)= [ I ■ F{dx) = 1. 
This completes the proof. □ 

4.2 Global Clustering Coefficient 

The global clustering coefficient is defined by 



1 " 

Cn — — / Cn{i)- 
n ^ — ^ 



z— 1 

Since it is a symmetric function of (.ti, . . . , a;„), we can prove SLLN for C„ by using the ergodic theory. 
Theorem 3. As n ^ 00, 

Cn E [C(l)] , almost surely. 

Proof. For simplicity, we only deal with the case E [C(l)] — 0. For general cases, we can prove the 
theorem by replacing C(f) by C(f) — E [C(l)]. Let x — {xi,X2, ■ . ■) be an infinite vector and Xk = Xk. 
We define measure-preserving transformation T„ for the product measure P by 

Xk+i if 1 < fc < n — 1, 
(T„a;)j, ^ {xi if fc = n, 
Xk otherwise, 

for each n > 1. By denoting C„(i; x) = C„(i; x^), a realization of C„ is represented by 

n 



n— 1 

C„(a;) = -VC„(l;r». 

n ^ — ^ 



i=0 



For arbitrary fixed e > 0, we define 

Cl(\;x)^{Cn[\-,x)-e)-lA^_. (10) 

ri-l 
1=0 

where \x : limsup^^o^ Cn{x) > e} . Using the maximal ergodic theorem (see Theorem 24.2 of [3]), 

Q(l;a;)(flP> 
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for every n> 1, where = {x : sup]^<j<„ Sj{x) > O}. On the other hand, we have 



S^(x) 

X : sup Sl{x) > \ = ix : sup ^ 



fe>i 



k>l 



> 0} = <x : sup C„(cc) > e} nA^ ^ As, 



k>l 



as n CX3 by Eq. ()10p . From the dominated convergence theorem and Theorem [2 we derive 

0< / Q(l;cc)(iP-> / [C{l;x)-e]dP, (11) 

as n ^ oo. 

Let Xn be the class of sets that are invariant under aU permutations of the first n coodinates and 
1 = rt^=i^n- It is easy to check that A^- £ X. Since P{A) equals to zero or one for any A e X by the 
Hewitt-Savage zero-one law (see Theorem 36.5 of [3]), the conditional expectation E[C(1)|X] equals to 
E [C(l)], almost surely. This leads to 



0< / [C{l;x)-e]dP = [ C{l-x)dF-eF{A,)= [ E[C(l;a;)|X] 
Ja, Ja, Ja^ 

E [C(l)] dV - eP{A,) ^ -eF{A,), 



dF-e 



by Eq. HIl) and E [C(l)] = 0. Then we have PiA^) = for any e > and therefore limsup„^oo C„ < 0, 
almost surely. Repeating the same argument for — C„, we have liminf„_^oo C„ > 0, almost surely. This 
completes the proof. □ 

Here we show a simple example for Theorem [3] Consider an i.i.d. sequence Xi, . . . ,X„ such that 
F{X, = 1) = p and F{X, = 0) = 1 - p for alH = 1, . . . , n. Let Bi = (6*, oo) and f^x, y) =x + y. We set 
a threshold 9 such that < 6* < 1 . In this case, a pair of vertices i and j with i 7^ j is disconnected if 
and only if Xi = Xj = 0. By direct computation, we have 



E[C(1)] = p-C(l;l) + (l-p)-C(l;0)=p 



+ 2p(l - p) 
1 



+ (l-p).i^ = l-p(l-p)^ 



In order to calculate C„, let Sn — J2^=i-^ij that is, the number of vertices with Xi = 1. We use the 
symbols Xi, s„ and c„ as realization of random variables Xi, Sn and C„ respectively. If Sn — 0, the graph 
consists of n isolated vertices. In this case c„ — w. If s„ — 1, the graph is the star in which only one 
central vertex has n — 1 edges and other n — 1 vertices are connected only to the center. So we obtain 



c„--{0-l + w(n-l)}= ( l-- 



If 2 < s„ < n — 2, s„ vertices with Xi — I have n — 1 edges, and the other n — s„ vertices are connected 
only to the vertices with Xi = I. So we have 



1 f i'V) - ("T") 



(n - Sn) > = 1 



(n - .S„)(n - 1 - Sn)Sr. 



If Sn = n — 1 OT n, we obtain the complete graph, and c„ = 1. Noting 

1 



{n~ l)(n- 2) 



1 



{n - Sn){n - 1 - Sn)Sn 
n(n- l)(n-2) 



1-i 



1 - 



{n — Sn){n — l — Sn)Sn 

ri(ri-l)(n-2) 



if Sn — 0,n — l,n, 
if Sn = 1, 
otherwise. 



we have 



Cn = 



[n - Sn){n - I ~ Sn)Sn 



1 - 



1-1- 



i{n- l){n-2) 

Sn 



{l{2,...,n}^Sn) + r ■ /{0,l}('5'n)} 



Sn 



1 - 



n - 1 



Sn 

n-2 



{l + (r-l)./{o,l}(^n)} 



1 - (1 -pYp = E [C(l)] , almost surely [n 00). 



8 



The last convergence conies from SLLN for the i.i.d. sequence. 

One of our motivations to study hmit theorems for the clustering coefficients is to make a clear 
distinction between the proportion of triangles in an entire graph and the clustering coefficients. By 
Eq. (66) of [9], the normalized number of triangles including vertex 1 converges to Et{x), almost surely 
for each realization x of Xi, where the normalization constant is equal to ("2^). In the same way, the 
degree of vertex 1 normalized by n — 1 converges to Ed{x), almost surely. Thus, when Ed{x) > 0, the 
local clustering coefficient converges almost surely to Et{x)/ Ed{x)'^ , that is, the limit of the normalized 
number of triangles divided by the square of the limit of the normalized degree. The mean field result 
corresponding to Theorem [2] is found in Eq. (3) of [18]. The denominator equals to Et{x) and the 
numerator [k{x)/N]'^ converges to Eoix^, almost surely as — * 00, where k{x) is the degree of the 
vertex 1 and N denotes the number of vertices. Equation (30) of [5] corresponds to the normalized 
number of triangles. These heuristic results are consistent with our rigorous result. Several examples for 
the global clustering coefficient are calculated in [10]. 

In practice, we may substitute into w and consider 

1 " 
n — number of vertices with degree or 1 '^"^^■^ 



i=l 



"-E.=i^{o}(K.(i)) ^ 
instead of C„. Using the same arguments of Theorems [5] and O it is easy to prove that 

^ n 1 

lim -y /.o}(K(z)) = lim -y/{^}(C„(i))=P(C(l) = u.)=P(£;i5(Xi)=0), almost surely 

1=1 i=l 

The last equality follows from the definition of C(l), i.e., Eq. ([3]). Noting that 

^"=l-(l/n)Er=i^{o}(K(*))'^"' 

we have the following: 
Corollary 2. As n ^ 00, 



5 Examples 



In this section, we show examples of the limit degree distribution, i.e., m — 2 and A2 is chosen as the 
collection of all possible edges in the limit theorem (Fact [3]) . We consider the case I = 1 with fc = fl- 
We assume that the random variable Xi is absolutely continuous so that it has a probability density 
function /. Let supp/ = {x G M : /(x) 7^ 0} be the support of /. 

We first set Bi = {9, 00) and fc{x, y) — x + y, i.e., the threshold network model in which an edge 
forms ii 9 < Xi + Xj for a given threshold 6* S R [5,6,9,10]. By calculating the characteristic function of 
D ~ U{Bi, A2), namely, the density of edges connected to vertex 1, we obtain the following results: 

Lemma 2. 1. If there exists a e M such that supp / = [a, 00), then 



5i{dk) if 9<2a, 

/(g-F-i(l^fc)) 

/(F-i(i-fc)) 

+ F[e - a)) ■ 5i{dk) if 9>2a. 



^ ~ S I{i-F{e-a)s)(k) ■ ,A dx 
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2. If there exists 6 e K such that supp/ = (—00, h], then 



' {l-F{0-b)) -Soidk) 

D ~ <; +/(o,i-F(.-.))(fc) ■ ^'['(^/_~[^_~J^^ dk if e<2b, 

^do{dk) if 0>2b. 

3. If there exist a,b gR such that supp / = [a, b], then 

'Siidk) if e<2a, 

li^zlll^dk 



D ~ i 



if 2a < < a + b, 
if = a + b, 



/(l-F(.-a),l)(fc)-^^^7T^ 

+ {l-F{0-a)) -Siidk) 

IiQ i)(k) ■ ^—7 v-^dk 

(l-F{0-b)) - Soidk) 

+J(o,i-F(^-6))(fc) • if a + b<0<2b, 

So{dk) if 6i>26. 

Furthermore, if f is symmetric on supp /, then 

£>~ J(o,i)(fc)rffc if e = a + b. 

4- If supp f = (—00,00), then 

n r /(^ - ^^"^1 - fc)) 

-D ~ /en 11 (fc) • — 7^ rr-^dk 

(0-1)^ ^ /(F-i(l-A;)) 

/or any ^ e K. 

Example 1. (Exponential distribution) If the random variable X\ has the probability density function 



fix) 



Xe-^^ if x>0, 
otherwise, 



(12) 



for a given A > 0, then 



D 



Siidk) 



if 0<O, 



I{e->-os)ik)-'^dk + e-^^ -diidk) if 9 > 0. 
Example 2. (Pareto distribution) If 



fix) 



'i-{^y^' if x>a, 
otherwise. 



for given a,c> 0, then 



D 



diidk) 



if e< 2a, 



hi^J^Mk)-{j:^Y^'dk+{^S.6,idk) if 0>2a. 



The distribution of D of these two examples is proportional to The exponent a equals 2 in 

Example 1 and 1 + 1/c in Example 2. Because of a lower cutoff of / in both examples, the limit 
distributions have weights on 61. 
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Example 3. (Uniform distribution) If 
then 



1 if 0<x<l, 
otherwise, 



6i{dk) if 0<O, 

I(i-e,i){k)dk + {l-e)-5i{dk) if < 9 < I, 

D ^ I I(^,^){k)dk if e = i, 

{9 - 1) ■ Soidk) + I(^o^2-e){k)dk if 1<9<2, 

So{dk) if 9>2. 

In this case, the limit distribution is mixture of the uniform distribution and the delta measure. 

By choosing Bi — {di, ^2] and /c(a;, y) = x + y, we obtain a generalization of the model investigated 
in [5,6,9,10]. More precisely, an edge {i,j) forms if 9i < Xi + Xj < 92 for given thresholds 6*1, G K 
such that 9i < 6*2. To calculate the characteristic function oi D ~ U{Bi, A2), we consider the case in 
which a random variable Xi has the probability density function (112p . i.e., the exponential distribution, 
for which the limit distribution is represented by: 



D 



Soidk) 



Soidk) 



+-^(0,l-e-^«2)(fc) 



(1-fc) 



-dk 



if 9i<e2< 0, 

if 6I1 < < 92, 



Soidk) 



"-^(o,i-e-^(»2-<'i))(fc) • (i-ky^ dk 



fe2 



-dk 



if < 6I1 < 



Finally, we deal with an example with 1 = 2. For fixed 
{e,(X)],B2 = (0,c], flix,y) = X + y, and flix,y) = \x ~ y 
distributed according to the exponential distribution (Eq. (fT^ ). This is the model proposed in [12] 
Because the kernel function of this model is 



e M and c £ [0, 00), we choose Bi = 
We consider the case in which Xi is 



/r 



— c+a:,c+ 



hix, X2) 



t{e-x,c+x] 



x]i^2) if '-^<x, 

X]{X2) 







if 

if a; < 



^ <- <- 6+c 
2 :i 2 ' 

e-c 
2 ' 



the limit distribution D 



UiCg^c,A2) is the following: 

-(l_e-M9-c)/2) ^^(^fc) 

+-^(0,2e-^(»+<^)/2 sinh(Ac)] (k) ' gik)dk 

rdk 



D 



if c< 



/(0,e-Ae_e-Ac](fc) • 2sinh(Ac)' 

+/(e-Ae_g-Ac i_g-A(8+<:)] (fc) • g{k)dk 

+ /(l,e-Me+c),l_e-^^'=](fc) ' 2 sin^ Ac) < < C, 



'(0,1- 



](^) ■ 2sinh(Ac) '^'^^^ 
:,l_e-2A<:](/c) • 2sinh(Ac) 



'(0,1- 



,dk 



'(1-e- 



2sinh(Ac 
,l„e-2Ac](/c) • 2sinh(Ac) 



if 



-c< 6* < 0, 



where 



4e- 



-AS 



if 6* < 



2 sinh(Ac) 
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